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含 单 裂隙 岩 体 蠕 变 破裂 裂纹 扩展 演化 规律 研究 


赵 娜 , 孟 利 新 , 张 怡 试 


(辽宁 工程 技术 大 学 力学 与 工程 学 院 ,123000 阜新 ) 


摘 要 : 岩 体 内 部 存在 大 量 不 同 尺寸 裂隙 ,形成 不 同 结构 形 式 的 裂隙 岩 体 ,裂隙 的 存在 会 降低 岩 体 
的 力学 性 能 ,对 于 岩 体 工程 的 整体 稳定 性 具有 很 大 的 影响 。 本 研究 基于 滑动 翼 裂 纹 模 型 对 恒定 载 
荷 下 含 单 裂隙 岩 体 蠕 变 变 形 进 行 了 理论 推导 和 分 析 , 基 于 MAPLE 软件 直观 给 出 恒定 载荷 下 裂隙 角 
度 、 长 度 对 里 变 变 形 的 影响 ,采用 连续 - 非 连续 数值 分 析 软 件 CDEM ,对 预制 单 裂 陈 在 不 同 倾角 、 不 同 
长 度 情况 下 岩 体 的 蠕 变 破裂 裂纹 扩展 进行 数值 模拟 。 通 过 理论 分 析 和 数值 模拟 可 以 看 出 , 岩 体 里 
变速 率 随 着 裂 阶 长 度 的 增 大 而 增 大 ; 随 着 裂隙 角度 的 增加 出 现 先 增 大 后 减 小 的 趋势 ,其 中 裂隙 倾角 
为 30" 时 蠕 变速 率 最 大 。 含 单 裂隙 岩 体 蠕 变 破 裂 裂 纹 扩 展 规律 的 研究 可 为 岩 体 工程 蠕 变 破裂 规律 
研究 及 蠕 变 治 理 提 供 一 定 的 理论 依据 和 参考 。 
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Study on the meso-crack propagation and evolution law of creep 
fracture in rock mass with single fissure 


ZHAO Na,MENG Lixin, ZHANG Yibin 


( College of Mechanics and Engineering, Liaoning Technical University , 123000 Fuxin , China) 


Abstract: There are multiple cracks of different sizes in rock mass, forming cracked rock masses of 
different structural forms. The existence of cracks not only reduces the mechanical properties of rock mass, 
but also significantly affects the overall stability of rock mass engineering. Based on the sliding wing crack 
model,this paper theoretically deduces and analyzes the creep deformation of rock mass with single 
fracture under constant load. Based on the MAPLE software , the impact of the crack angle and length on 
the creep deformation under constant load is intuitively given. Apart from that , continuous-discontinuous 
numerical analysis software GDEM is used to numerically simulate the creep rupture crack propagation of 
rock mass under different inclination angles and different lengths of prefabricated single fissures. By 
theoretical analysis and numerical simulation, it can be observed that the creep rate of rock mass increases 
with the increase of crack length. Moreover, with the increase of crack angle,the creep rate first increases 


and then decreases, and it is the largest when the crack inclination angle is 30?. In short, the study on 
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creep rupture crack propagation law of rock mass with single fissure can provide certain theoretical basis 


and reference for research on creep rupture law and creep control of rock mass engineering. 


Key words :single crack ; creep rupture; crack propagation ; evolution law ; constant load 


岩石 变形 过 程 中 与 时 间 相 关 的 现象 称 为 流 变 。 
流 变 包括 蚂 变 、 松 弛 和 弹性 后 效 。 岩 石 里 变 是 指 在 
应 力 不 变 的 情况 下 ,变形 (或 应 变 ) 随时 间 增 长 的 现 
象 。 应 力 不 变 ,岩石 变形 随时 间 逐 渐 增 长 的 主要 原 
因 是 岩石 的 非 均 匀 性 。 岩 石 的 非 均 匀 性 包括 材料 非 
均匀 性 和 结构 非 均 匀 性 ,材料 非 均 匀 指 岩石 由 不 同 
的 矿物 质 元 素 组 成 ,而 结构 非 均 匀 性 指 岩 石 中 含有 
的 孔隙 、 裂 阶 和 孔洞 等 。 天 然 岩 体内 部 存在 着 大 量 
的 不 同 尺 寸 的 裂隙, 形成 不 同 结构 形式 的 裂隙 岩 体 。 
在 岩 体 工程 中 , 裂 际 孔隙 的 结构 形式 对 于 岩 体 工程 
的 整体 稳定 性 具有 很 大 的 影响 。 不 同 工 程 环境 下 的 
岩 体 结构 裂纹 扩展 形式 .变形 响应 及 规律 也 不 尽 相 
[n] ,这 很 大 程度 地 影响 了 岩 体 结构 工程 的 稳定 性 。 
含 裂隙 宕 体 的 裂纹 扩展 过 程 更 是 十 分 复杂 ,裂隙 的 
存在 会 降低 岩 体 的 力学 性 能 , 岩 体 内 部 裂隙 扩展 贯 
通 形成 宏观 破坏 面 。 由 此 所 引发 的 工程 失 稳 破坏 事 
件 屡见不鲜 。 影 响 裂 孙 岩 体力 学 性 能 的 因素 很 多 ， 
例如 : 岩 体 裂隙 的 尺寸 .间距 ,倾斜 角度 等 。 国 内 外 
专家 和 学 者 对 此 已 做 过 一 些 研 究 。 杨 圣 奇 等 ' ”| 通 
过 单 轴 压 缩 试验 ,探究 了 含 裂隙 砂 岩 的 力学 特性 ,分 
析 了 裂 陀 倾 角 对 含 缺陷 砂岩 试 件 强度 和 变形 行为 的 
影响 。LI 等 中 通过 对 单 轴 压 缩 下 含 裂隙 的 大 理 岩 
试 件 裂 纹 扩展 进行 试验 ,结果 表明 含 预制 裂隙 试 件 
的 临界 失 稳 载荷 、 破 坏 时 的 应 变 以 及 弹性 模 量 都 有 
所 降低 。NOLEN-HOEKSEMA 等 采用 对 真实 岩石 
人 为 预制 裂隙 的 方法 研究 裂纹 对 岩石 的 影响 ,在 大 
理 岩 上 预制 斜 裂纹 并 进行 了 单 轴 压 缩 试验 ,人 研究 发 
现 裂隙 尖端 产生 的 新 裂纹 不 仅 有 所 扩展 而 且 扩 展 形 
式 不 以 裂纹 为 对 称 轴 对 称 。 此 外 ,还 有 更 多 学 者 在 
基于 单 裂隙 的 基础 上 对 多 裂隙 岩 体 进行 了 研究 ,如 
机 东 伟 等 年 对 断 续 三 裂隙 岩 体 开展 单 轴 压 缩 实验 ， 
得 到 不 同 中 心 裂隙 倾角 下 的 涯 样 损伤 劣化 模式 。 薄 
成 志 等 ““" 采用 伺服 控制 的 单 轴 加 载 系 统 对 预制 多 
裂隙 的 水 泥 砂浆 试 件 进行 试验 ,探讨 了 裂隙 角 度 和 
分 布 密度 对 类 岩石 试 件 的 影响 ,发 现 不同 裂 隙 倾角 
对 试 件 峰值 强度 的 影响 不 同 。 相 比 单 轴 压 缩 试 验 之 
外 , 王 辉 等 对 含 袭 隙 试 件 开展 了 巴西 圆 盘 试 件 的 
臂 裂 试验 ,得 出 了 层 理 与 预制 裂隙 共同 作用 下 岩 样 
的 断裂 特性 和 破坏 机 制 。 赵 程 等 "利用 数字 图 像 


相关 方法 (DICM) ,对 含 不 同 预 置 单 裂 纹 的 类 岩石 材 
料 在 单 轴 压 缩 作 用 下 的 裂纹 扩展 规律 及 岩 体 细 观 损 
伤 演化 机 制 进行 了 系统 研究 ,研究 表明 随 着 裂隙 角 
度 的 变化 ,裂纹 形式 和 最 大 应 变 集 中 区 域 不 同 。 李 
树 刚 等 2 利用 自主 研发 的 YYW-I 应 变 控制 式 无 侧 
限 压 缩 仪 及 高 精度 声 发 射 (AE) 监测 系统 ,研究 了 
含 不 同 宽度 预制 裂隙 的 类 岩石 材料 试 样 在 单 轴 压 缩 
条 件 下 的 受 力 变形 和 破坏 模式 ,分 析 了 预制 裂隙 宽 
度 对 类 岩石 材料 破坏 的 影响 规律 。 然 而 ,与 实验 相 
比 ,数值 模拟 具有 成 本 低廉 ,处 理 和 分 析 结 果 快 捷 等 
优点 ,能 针对 同一 个 岩 样 进行 多 次 预制 变量 的 对 比 
试验 ,从 而 避免 了 实验 中 由 于 岩 样 不 同 造成 的 误差 ， 
从 而 保证 试验 结果 的 准确 性 。 张 仕 林 等 “采用 单 
轴 压 缩 试验 ,结合 数值 模拟 方法 ,研究 了 红 砂 岩 裂 际 
在 不 同 张 开 度 条 件 下 的 扩展 机 理 。XU A5 P 
轴 和 双 轴 压缩 数值 试验 ,发 现 岩 石 破 坏 演化 首先 发 
生 裂隙 周围 区 域 。HUANG 等 "采用 数值 模拟 方式 
对 含 双 裂隙 砂岩 实验 结果 进行 裂纹 起 牧 位 置 及 扩展 
方式 的 验证 。 

上 述 成 果 丰 富 了 含 裂隙 岩 体 变形 破坏 方面 的 研 
究 , 对 含 裂隙 岩 体 袭 纹 的 产生 以 及 含 预制 裂隙 与 应 
力 之 间 的 关系 进行 了 研究 ,而 关于 在 恒定 外 载荷 作 
用 下 裂纹 随时 间 二 次 扩展 直至 贯通 发 生 宏观 蠕 变 破 
裂 的 研究 仍 比 较 薄 弱 。 本 研究 通过 对 恒定 载荷 下 合 
单 裂隙 岩 体 的 蠕 变 变 形 理论 推导 ,给 出 恒定 荷载 下 
裂隙 角度 .长 度 对 蠕 变 变形 的 影响 ,采用 连续 - 非 连 
续 数 值 分 析 软 件 CDEM , 对 预制 单 裂隙 在 不 同 倾角 、 
不 同 长 度 情 况 下 岩 体 的 宏 细 观 破裂 机 理 进 行 研 究 ， 
为 裂隙 岩 体 蠕 变 破裂 裂纹 扩展 规律 研究 及 蠕 变 治理 
提供 理论 依据 。 


NC 


1 单 裂隙 岩石 蠕 变 变形 破裂 


含 初始 单 裂 际 宕 石 在 单 轴 压缩 状态 下 的 受 力 简 
化 模型 如 图 1(a) 所 示 。 根 据 线 弹 性 断裂 力学 可 知 ， 
当 裂 际 尖 端的 应 力 小 于 临界 应 力 (断裂 韧 度 ) T, A 
纹 尖 端 不 会 发 生 裂 纹 扩展 ; 当 尖 端 应 力 大 于 或 者 等 
于 临界 应 力 ( 断 裂 韧 度 ) AF, RRR mT E RNA 
纹 , 在 两 个 初始 裂纹 端 部 发 生 裂纹 扩展 。 
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图 1 单 裂 除 岩 石 简化 模型 
Fig.1 Simplified model of single fractured rock 
AUSUR BUE nd 5€ HARE SUE TE rp h BRUST TREE 
坏 的 一 种 裂纹 形式 , d rh ERAD HE TRE JT VETT TEE o 


TRUST JE PT TEE D e S Ee fa 
d:ptbepim sacas m. odd 
WAHIEROA IRRA, 裂纹 长 度 为 254 村 
AAT Ke BER RADIAN o T A CURE 
系数 为 BRADH Fusu 为 常数 。ASHBY ^87 4 
H AEREE RR A R LA INTR. 
裂纹 尖端 的 工 型 应 力 强 度 因子 可 表示 为 
Kuce cl (1) 
Ln(1 Bb) 1* 
tH e e RRAN HS ERE n 
向 一 致 。 作 用 于 初始 微 裂纹 上 的 独 人 力 FI e 
维 条 件 下 的 表达 式 可 写 为 


F, = [e + ug) Tb sina --A,ob (2) 


式 中 ,4 为 与 摩擦 系数 B ORC u ARSA, EU 


4 = m JEY +w +u) (3) 
由 断裂 力学 可 知 , 在 不 考虑 亚 临 界 裂纹 扩展 作 
用 的 情况 下 ,当主 裂纹 产生 实际 滑动 , 即 作 用 于 主 裂 
纹 上 的 切 应 力克 服 摩擦 阻力 时 ,其 裂纹 尖端 才 进 行 
扩展 。 此 时 临界 应 力 强度 因子 记 为 Kic ,临界 应 力 可 
E. 27:9 
4 Kic 


Tic =- EE (4) 
i: 14g -u mb 


鉴于 亚 临 界 裂纹 扩展 的 影响 ,裂纹 在 Ki < Ki 
时 ,脆性 岩石 裂纹 开始 生长 。 对 于 微 裂 纹 ,裂纹 增长 
Ax] gogo m 


zh E n. 为 应 力 腐蚀 指数 ;> 为 特征 裂纹 速度 。 


H 力学 学 报 


第 40 卷 
两 个 稼 数 都 可 以 通过 亚 临 界 裂纹 扩展 测试 获得 。 将 
式 (1) 代 入 式 (5) ,可 以 得 到 随时 间 变 化 的 必 有 裂纹 长 
度 的 微分 方程 为 

di - v (sb)? [Arke a (6) 


dt 
Xm, 


IC 


e, = m? (I/b +B)” (T) 
& 22 (mb)? Vr [ D^ - (1. LJ] (8) 


脆性 岩石 中 往往 存在 许多 初始 微 裂纹 ,并 且 微 
裂纹 的 分 布 是 随机 的 。 根 据 WeiBull 概率 分 布 '2) ， 
岩石 损伤 可 表示 为 

€ 

D-i-ex|-(£) | (9) 
IP: m 和 so 是 材料 常数 ; es 是 轴 向 应 变 。 由 
BRISTOW” 定义 且 被 WALSH 9 推广 的 关于 初始 
弄 纹 密度 对 岩石 的 力学 行为 可 知 , 含 单列 隙 岩石 初 
始 损 伤 可 定义 为 

D, = 4TN (ba) /3 (10) 
式 中 ,N 是 初始 单 裂纹 的 单位 体积 值 ,将 初始 损伤 代 
人 方程 式 (9) ,可 以 得 出 裂纹 扩展 之 前 的 弹性 应 


变 为 


e, = e [-Inc1 - Do1"7" (11) 
随 着 裂纹 的 扩展 , 微 裂 纹 的 密度 和 大 小 发 生变 
化 ,损伤 可 以 定义 为 ”| 


D = AN, rab (12) 


比较 式 (10) 与 式 (11) , RE DAS HUZDUL I ZEE 55 3E 
裂纹 的 成 长 之 间 的 关系 为 
e =e0[-ln[l- (+ (CoD Di (Cab 
(13) 
里 变 失 效 随时 间 变 化 的 经 裂纹 长 度 ( 即 100) ) 
可 以 为 由 式 (6) 的 数值 积分 进行 计算 。 用 演化 解 代 
替 辟 裂纹 长 度 代 人 式 (13) ,可 以 得 出 应 变 和 时 间 之 
间 的 关系 为 
e (t) 2e [ -In [1 - (L(t) +(ab)) Do/(ab) ] ]" 
(14) 
式 (14) 可 以 反映 某 一 岩石 试 样 ,相同 的 破坏 值 微观 
尺度 和 宏观 尺度 之 间 的 关系 。 
3e D LII) (L(t) 十 Cabil- 
m[1 - (I(t) + Cab)’ Dl 


elt) = 


[- i [i - 9 e canna” 
(15) 


投稿 网 站 :http://ejam. xjtu. edu. cn” 微 信 公众 号 :应 用 力学 学 报 


第 5 Hj 


对 式 (14) 进 行 时 间 的 求 导 可 得 式 (15 ) ,得 到 岩 
石 试 样 里 变速 率 与 变量 之 间 的 关系 。 

为 了 反映 微观 裂纹 矿 才 与 宏观 应 变 的 关系 , 基 
于 滑动 辟 刚 纹 模 型 ,假设 在 给 定 的 临界 微 裂 纹 长 度 
下 会 发 生 宏观 破坏 ,采用 maple 软件 可 得 出 式 (14) 
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中 不 同 裂 际 倾角 和 不 同 裂 际 长 度 情况 下 的 岩 体 应 变 
随时 间 的 关系 ,但 基于 滑动 疲 裂 纹 模型 研究 假设 的 
局 限 性 ,maple 软件 只 能 表达 出 岩 体 蠕 变 的 初始 蠕 变 
阶段 和 稳定 蠕 变 阶 段 , 如 图 2 ~ 3 所 示 , 给 定 参数 如 
表 1 ,参数 均 取 自 于 中 国 锦 屏 大 理 石 的 实验 结果 。 


表 1 参数 取 值 
Tab.1 Parameter values 
Kie/ ( MPa ml) lmm (m.s?!) D, B £o " " n 
1.61 2 0.16 0.048 0.32 0.014 7 57 0.51 1 
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Fig.2 Creep curves of fractured rock mass with 
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Fig.3 Creep curves of fractured rock mass at different angles 


图 2 和 图 4 nE], b 25 mm 时 , 岩 体 的 初 
始 蠕 变 阶 段 总 应 变 为 1. 42 x 10 7 fa a Wa D Ex 
变速 率 为 1.24 x10 5/5; 4 b 26 mm 时 ,初始 蠕 变 阶 
段 总 应 变 为 1. 52 x 10 7 , 稳 态 蠕 变速 率 为 1. 283 x 
10 5/s;?4 b 27 mm 时 ,初始 蠕 变 阶 段 总 应 变 为 
1.56 x 107 ,稳定 里 变速 率 为 1.32 x 10 "7s; HE 
可 见 , 含 不 同 初 始 长 度 和 裂隙 岩 体 蠕 变 曲 线 趋势 大 致 
相同 , 岩 体 的 蠕 变速 度 随 着 裂隙 的 长 度 变 化 而 发 生 
相应 变化 , 裂 际 长 度 越 长 ,里 变速 率 越 快 。 

由 图 3 和 图 5 可 知 , 当 a=30° 时 ,里 变 初始 蠕 变 
阶段 总 应 变 为 5.46 x 10 ,稳定 蠕 变 阶 段 蠕 变速 率 


为 1.36 x 10 7/5; 5 a =45° 时 ,初始 蠕 变 阶段 总 应 
变 为 1.42 x10 ,稳定 蠕 变速 率 为 1.24 x 10 575; 14 
a =60° 时 ,初始 蠕 变 阶 段 总 应 变 为 0.62 x10 ,稳定 
里 变速 率 为 1. 12 x 1075/5; 24 a =90* 时 ,初始 蠕 变 
阶段 总 应 变 为 0. 22 x 10 7 , 稳定 蠕 变 速率 为 0.96 x 
10“”/s。 可 见 含 单 裂隙 岩 体 在 不 同 裂隙 倾角 的 情况 下 
岩 体 蜂 变 速率 也 不 尽 相 同 ,倾角 过 大 过 小 都 会 影响 兰 
体 蠕 变速 率 , 当 倾角 为 30° 时 ,里 变速 率 叶 现 出 最 快 。 


1.58 
1.56 
1.54 初始 阶段 总 应 变 
30 . 
1.52 Lo 
ja x 
S 1.50 2 
Bi 148 pa 
18 x 
1.46 dh 
1.44 
142 
140 
5.0 5.5 6.0 6.5 7.0 


长 度 /mm 
图 4 不 同 长 度 裂 际 岩 体 初始 阶段 总 应 变 与 稳 态 里 变速 率 
Fig.4 Total strain and steady-state creeprate in initial stage of 


fractured rock mass with different lengths 


EIE 


总 应 变 /10” 
稳 态 里 变速 率 /10” 


30 45 60 90 
HEC) 
图 5 不 同 角度 裂隙 岩 体 初始 阶段 总 应 变 与 稳 态 


Fig.5 Total strain and steady state in initial stage of fractured 


rock mass with different angles 
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2 ”里 变 破裂 裂纹 扩展 演化 数值 模拟 


2.1 连续 - 非 连续 变形 计算 基本 原理 


连续 - 非 连续 变形 计算 是 基于 连续 介质 力学 的 
离散 元 方法 ,该 方法 将 有 限 元 与 离散 元 进行 耦合 ,在 
块 体内 部 进行 有 限 元 计算 ,在 块 体 边界 进行 离散 元 
计算 ,不 仅 可 以 模拟 材料 在 连续 状态 下 及 非 连续 状 
态 下 的 变形 .运动 特性 ,更 可 以 实现 材料 由 连续 体 到 
非 连续 体 的 渐进 破坏 过 程 。 

连续 - 非 连续 模型 包括 由 微 元 素 组 成 的 弹性 耦 
合 表 征 元 素 。 由 于 其 组 成 微 单元 的 断裂 条 件 不 同 ， 
每 个 表征 单元 具有 不 同 的 弹性 模 量 (图 6) 。 对 于 给 
定 的 表征 元 素 ,弹性 模 量 分 布 是 根据 其 内 部 的 断裂 
程度 分 布 计算 得 出 的 ,反映 了 岩石 的 微观 非 均 质 性 。 
表征 元 素 之 间 的 连接 (弹性 或 断裂 ) 模拟 了 微观 裂 
颖 的 萌生 .扩展 和 连接 ,因此 代表 了 在 渐进 式 岩 石 蠕 
变 和 破裂 过 程 中 宏观 裂缝 的 演变 


表征 元 


弹性 微 元 ” c 
弹性 微 元 破裂 微 元 


图 6 连续 - 非 连续 计算 原理 


Multiscale rock model and its relations 


Fig. 6 


2.2 计算 模型 及 方案 


根据 理论 分 析 可 以 看 出 ,在 不 同 的 裂隙 长 度 、 裂 
隙 倾角 下 , 岩 体 蠕 变 变 形 特性 不 同 。 为 更 好 的 研究 
岩 体 蠕 变 过 程 中 裂纹 扩展 演化 规律 ,采用 连续 - 非 连 
续 软件 CDEM 对 含 单 裂 际 岩 石 在 恒定 载荷 作用 下 
的 裂纹 扩展 过 程 进行 数值 模拟 ,分 析 在 不 同 的 裂 际 
长 度 、 裂 隙 倾角 下 岩 体 的 蠕 变 曲 线 以 及 裂纹 扩展 演 
化 规律 。 计 算 模型 如 图 7 所 示 ,a 为 裂隙 与 水 平方 
向 的 倾角 ,b 为 裂 际 的 半 长 度 。 模 型 尺寸 为 50 mm x 
100 mm ,选取 不 同 裂 际 角 度 和 长 度 进行 计算 ,计算 方 
案 如 表 2。 
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图 7 数值 计算 模型 
Fig.7 Simplified model diagram of test piece 
R2 计算 方案 


Tab.2 Calculation scheme 


裂 际 角度 A/(°) 裂隙 长 度 /mm 
方案 一 0 10 
15 10 
30 10 
45 10 
60 10 
90 10 
方案 二 45 10 
45 12 
45 14 
连续 - 非 连续 变形 破裂 过 程 中 ,连续 变形 和 单元 


内 采用 Mohr-Coulomb 模型 ,参数 设置 如 表 3。 
表 3 单元 内 计算 力学 参数 


Tab.3 Calculated mechanical parameters of rock mass 


材料 密度 / 抗 拉 强 度 / 内 摩擦 角 / 弹性 模 量 / 粘 聚 力 /” 剪 胀 角 / 泊 松 
(kg*m^?) MPa (9) MPa MPa (?) 比 
2 500 3 40.0 1 x104 10 10.0 0.25 


对 于 连续 - 非 连续 模型 块 体 离散 元 ,单元 和 单元 


之 间 采 用 可 变形 表征 元 接触 弹簧 的 虚 接 触 面 , 弹 簧 
连接 的 本 构 模 型 为 脆 断 模型 (britleMC ) ,材料 参数 
如 表 4 所 示 。 

表 4 岩 体 接触 面 力 学 参数 


Tab.4 Mechanical parameters of rock mass contact surface 


法 向 刚度 / ” 切 向 刚度 / ARJ 抗 拉 强 度 / ”摩擦 角 / 
(Pa * m) (Pa * m) MPa MPa (S) 

5 x 10? 5 x 10? 10 0.4 30.0 
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2.3 岩石 蠕 变 破裂 裂纹 扩展 演化 过 程 


建立 单 裂隙 岩 体 简化 模型 , 裂 际 长 度 为 10 mm, 
倾角 为 30° ,采用 连续 - 非 连续 计算 软件 GDEM 建立 
模型 并 进行 网 格 划 分 ,网 格 单元 为 三 节点 的 三 角形 
单元 ,模型 底部 固定 ,根据 实验 经 验 对 顶部 施加 恒定 
应 力 4 MPa ,数值 模拟 时 间 步 长 为 dd =1 s。 为 了 获 
得 岩石 恒定 压力 下 的 蠕 变 曲线 ,在 模型 上 边缘 中 点 
设置 位 移 监 测 点 ,根据 监测 点 数据 绘制 蠕 变 曲线 。 
蠕 变 曲 线 图 如 图 8 所 示 , 不 同 蠕 变 阶段 裂纹 扩展 演 
化 过 程 如 图 9 所 示 。 

从 图 8 和 图 9 可 以 看 出 , 含 单 裂隙 岩 体 在 恒定 
REMEH F , 细 观 裂纹 萌生 .扩展 ,宏观 表现 为 蠕 变 
变形 。 岩 石 蠕 变 共 分 为 3 个 阶段 ,减速 蠕 变 阶 段 , 岩 
石 发 生 弹性 变形 , 随 着 时 间 的 增加 ,裂纹 在 初始 裂隙 
两 端 产生 ;等 速 蠕 变 阶 段 ,裂纹 进一步 扩展 、 延 伸 , 宏 
观 变形 增 大 ;加 速 蠕 变 阶 段 , 腊 纹 贯通 ,形成 宏观 断 
裂 面 ,此 时 岩石 仍然 具有 一 定 的 承载 能 力 , 此 后 宏观 
变形 加 速 , 蠕 变速 率 逐 渐 增 大 。 

从 图 9 裂纹 扩展 过 程 可 以 看 出 ,在 0 步 到 1200 
步 过 程 中 岩 体 处 于 弹性 阶段 ,无 次 生 裂 际 产生 ; 


0 步 1 200 步 1 250 步 2 900 步 


1 
弹性 阶段 i 减速 蠕 变 阶 段 


赵 娜 ， 等 : 含 单 裂隙 岩 体 蠕 变 破裂 裂纹 扩展 演化 规律 研究 


5 550 步 


1111 


1 200 步 初始 辟 裂 纹 a 在 预制 裂 际 的 下 端 产生 ,并 垂 
直 于 预制 裂隙 ;在 1 250 PRR b 在 预制 裂 险 上 端 
产生 ;至 2 900 步 时 裂纹 a 持续 扩展 为 裂纹 de; 当 
到 5 550 PRR Fm ERER c; 在 6 400 步 裂 
纹 持 续 扩 展 产 生 新 的 裂纹 g f; 当 扩展 到 10 000 步 
时 ,预制 裂纹 上 端 开 始 产生 次 生 裂 纹 ,并 在 40 000 步 
与 裂纹 h 相 联通 形 成 新 的 裂纹 j ,在 60 000 步 时 ,最 
终 裂纹 贯通 ,岩石 试 件 发 生 明 显 的 宏观 破坏 带 , 辟 裂 
纹 的 扩展 方向 与 恒定 载荷 方向 一 致 。 
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减速 ! 等 速 加 速 
0.000 01H KEZE: Mire 蠕 变 
阶段 阶段 WE 
0 — 19 00020 00030 00040 000, 00050 00070 00080 000 
时 步 
图 8 含 单 裂 隙 宕 体 蠕 变 变形 曲线 


Fig.8 Creep curve of rock mass with single fracture 
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Fig.9 Creep rupture crack propagation process of rock mass with single fissure 


2.3.1. RIRA RITE A I KARAI D 

选取 含 单一 裂 际 长 度 为 10 mm ,宽度 为 2 mm 的 
岩 体 模型 ,测试 了 不 同 裂 际 角度 (0°、15° 、30° 45°, 
60°* 和 90°) 下 的 单 轴 压 缩 蠕 变 实验 数值 模拟 结 
可 以 得 到 在 相同 材料 参数 及 加 载 应 力 情况 下 ,和 裂 际 
角度 变化 时 不 同 岩 体 蠕 变 破裂 及 等 效 位 移 云图 ,如 
图 10 所 示 。 

从 图 10 可 以 看 出 ,在 恒定 载荷 下 , 蠕 变 等 效 位 
移 和 发 生 蠕 变 断 裂 时 的 宏观 行为 对 于 不 同 的 断裂 角 
度 情况 有 显著 差异 。 当 a =0° 时 ,主要 变形 方式 为 
拉 剪 ,产生 垂直 恤 形 裂纹 ,倾斜 的 次 生 裂 纹 不 断 扩 大 
延伸 ,但 结构 相对 稳定 ,破裂 程度 较 小 。 当 a = 15? 


时 , 裂 际 尖 端 出 现 辟 裂纹 和 次 生 裂纹 , 岩 体 发 生 剪 切 
破坏 ,破裂 程度 较 a =0° 增 强 。 在 a =30° 的 情况 下 ， 
裂纹 尖端 同时 产生 到 形 裂纹 和 共 面 .倾斜 的 二 次 裂 
纹 。 裂 颖 不 断 发 展 , 经 过 一 定 程度 的 膨胀 后 ,通过 
“预制 "断裂 发 生 X 型 共 斩 剪 切 破 坏 。 对 于 a —45*, 
裂 颖 从 主 裂 颖 的 尖端 开始 , 辟 裂 颖 继续 扩大 , 直到 发 
生 穿 透 破裂 ,但 在 该 角度 ,预制 裂 颖 周围 的 破裂 效果 
比 a =30" 时 略 弱 。 在 a =60" 的 情况 下 ,裂纹 扩展 效 
应 似乎 减弱 , 沿 主 断 口上 端 纵向 渗透 。 当 a = 90? 
时 ， 裂隙 岩 体 形成 相对 稳定 的 结构 ， 在 荷载 作用 
下 ， 最 大 主 应 力 方 向 平行 于 裂 际 纵向 ， 裂 纹 扩展 不 
明显 。 
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图 10 不 同和 裂隙 角度 
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度 岩 体 及 等 效 位 移 云 图 


Fig. 10 Rock masses with different fracture angles and equivalent displacement cloud diagrams 


从 图 11 和 图 12 可 以 看 出 ,在 1200 25 zz Bij , PB 
裂 岩 体 的 特征 是 所 有 断裂 角度 的 瞬时 弹性 变形 。 从 
1 200 步 到 20 000 步 ,裂纹 萌生 发 生 在 不 同 的 方向 ， 
但 每 个 断裂 角 的 蠕 变 曲线 相似 ,为 减速 蠕 变 阶 段 。 
从 20 000 步 开始 ,变形 进入 等 速 蠕 变 阶 段 。 从 图 12 
可 以 看 出 ,此 时 不 同 断 裂 角度 的 里 变速 率 明显 不 同 。 
对 于 a =0。, 蠕 变速 率 为 6.05 x107/s; 对 于 a = 
15° 蠕 变速 率 为 8. 36 x 10 "7s; 对 于 a 2309, 蠕 变 
速率 为 1. 06 x 10 5/5; 对 于 a = 45? , IE AE xli RON 
7.03 x 10 77s; 对 于 a = 60°, 蠕 变速 率 为 6. 00 x 
10 ?/s;X[-F a =90。, 蠕 变速 率 为 4.87 x10 77s, K 
约 50 000 步 后 ,a =30° 的 情况 开始 出 现 加 速 蠕 变 ， 
进入 第 三 阶段 (加 速 蠕 变 ) 。 当 加 载 达到 约 60 000 
步 时 ,a =45",60" 上 断裂 的 情况 也 开始 进入 加 速 蠕 变 
阶段 。 总 体 而 言 , 岩 体 蠕 变速 率 随 初始 断裂 角 的 增 
大 先 增 大 后 减 小 。 断 裂 倾角 为 30" 时 蠕 变 率 最 高 ,在 


这 种 情况 下 , 岩 体 结构 容易 破裂 。 同 时 与 图 3 对 比 
分 析 可 得 到 理论 部 分 与 数值 模拟 中 相同 部 分 裂隙 角 
度 的 改变 ,对 含 裂 际 岩 体 的 蠕 变 破裂 曲线 影响 规律 
相同 。 
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图 11 不 同 角 度 裂隙 岩 体 宏观 蠕 变 曲线 图 


Fig. 11 Macroscopic creep curves of fractured rocks 
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稳定 里 变速 率 /(10”。 


角度 /(°) 


图 12 不 同 角 度 裂 际 岩 体 稳定 蠕 变 
Fig.12 Stable creep rate of fractured 
masses at different angles 
2.3.2 数值 模拟 与 实验 对 比 验 证 
连续 - 非 连续 软件 CDEM 用 于 模拟 在 恒定 载荷 

下 具有 单个 预先 存在 的 裂隙 的 岩石 中 的 裂纹 扩展 。 
在 材料 参数 和 加 载 应 力 相 同 的 情况 下 ,可 以 获得 岩 
体 模型 和 单 轴 压 缩 数 值 实验 。 不 同 的 初始 破裂 角度 
导致 不 同 的 岩 体 破裂 模式 ,这 些 可 以 与 蒋 明镜 等 ”” 
的 实验 室 结果 进行 比较 ,如 图 13 所 示 。 

室内 实验 和 GDEM 数值 模拟 均 表明 ,对 于 不 同 
的 初始 破裂 角度 a( 其 中 a 表示 裂纹 与 水 平面 的 夹 
角 ) 存 在 以 下 现象 。 

当 a=0° 时 ,站 裂纹 沿 预制 裂纹 尖端 沿 轴 向 载荷 
方向 扩展 , 随 着 时 间 的 推移 会 出 现 较 小 的 次 生 裂 纹 。 

当 aq=30° 时 ,裂纹 萌生 并 沿 预制 裂纹 扩展 ,最 
终 渗 透 形成 拉 伸 断 裂 。 右 豆 裂 纹 的 下 端 向 右 延 伸 ， 

贯穿 的 二 次 裂纹 。 

当 a=45° 时 , 主 裂缝 右 四 裂纹 扩展 并 穿 透 。 预 
制 裂 颖 的 右 下 端 与 次 生 裂 颖 形成 闭合 三 角形 面 。 
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14 不 同 裂 际 长 度 岩 体 及 等 效 位 移 云图 


Fig.14 Rock masses with different fracture lengths and equivalent displacement cloud diagrams 
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当 a =60° 时 , 剪 切 裂缝 在 预制 倾斜 裂缝 的 右 端 
萌生 ,并 沿 加 载 方向 扩展 直至 穿 透 。 


a 
a=0° 


a=308 a=45° a=60° 
图 13 GDEM 模拟 与 室内 试验 试 样 破坏 形态 对 比 图 
Fig. 13 Comparison of failure morphology of GDEM simulation 


and laboratory test specimens 

从 这 个 比较 中 可 以 清楚 地 看 出 ,裂纹 扩展 模式 
的 GDEM 模拟 与 试验 示例 非常 相似 ,其 优点 是 可 以 
更 容易 地 在 数值 上 跟踪 破裂 起 始 、. 扩 展 和 发 展 方向 。 
基于 此 ,连续 - 非 连续 变形 软件 CDEM 可 以 有 效 地 模 
拟 岩 石 裂 颖 扩展 的 演化 过 程 。 
2.3.3 有 裂隙 长 度 对 岩石 蠕 变 破 裂 的 影响 

对 具有 2 mm 宽度 的 单个 45 裂隙 的 岩 体 模型 
进行 了 数值 模拟 ,用 于 单 轴 压 缩 蠕 变 实验 ,其 中 测试 
了 不 同 裂 队长 度 (10 .12 14 mm) 的 恒定 材料 参数 和 
加 载 应 力 。 不 同 长 度 的 单个 裂纹 的 里 变 破裂 及 等 效 
位 移 云 图 如 图 14 所 示 ,不 同 长 度 的 单个 裂纹 的 里 变 
变化 曲线 如 图 15 所 示 。 
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图 14 表明 , 随 着 裂 际 长 度 的 增加 , 岩 体 空隙 体 
积 增 大 ,结构 不 均匀 性 增加 , 妥 纹 扩展 和 破裂 变 得 更 
加 严重 。 蠕 变速 率 与 裂 陈 长 度 呈 正 相关 。 
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图 15 不 同 长 度 单 裂 际 蠕 变 变 化 曲线 图 
Fig.15 Creep curve of a single crack with different lengths 
由 图 15 可 知 ,在 蠕 变 开始 时 ,不 同和 裂隙 长 度 的 
蠕 变速 率 是 相似 的 。 当 变形 进入 等 速 阶 段 时 , 蜂 变 
速率 受 裂 纹 长 度 的 影响 更 为 明显 :裂纹 越 长 ,里 变速 
率 越 高 ,进入 加 速 蠕 变 阶 段 越 时 。 蠕 变速 率 和 破裂 
程度 都 随 着 裂纹 长 度 的 增加 而 增加 。 通 过 与 图 2 不 
同 初始 长 度 裂 际 岩 体 蠕 变 曲线 图 对 比 ,可 得 到 在 理 
论 计算 和 数值 模拟 中 ,裂隙 长 度 的 改变 ,对 含 裂隙 宕 
体 的 蠕 变 破裂 曲线 影响 规律 相同 。 


2.4 不 同 变量 对 蠕 变 破裂 裂纹 扩展 影响 机 理 


对 于 含 单 裂隙 岩 体 蠕 变 破 裂 裂 纹 的 扩展 在 裂 际 
两 个 端点 处 萌生 并 进一步 扩展 ,在 宏观 表现 为 蠕 变 
变形 。 裂 纹 的 扩展 演化 与 裂 际 的 角度 和 长 度 密切 相 
关 。 通 过 以 上 含 单 一 裂隙 岩 体 在 不 同 变量 情况 下 的 
数值 模拟 分 析 , 可 以 得 出 对 于 裂 际 角度 (a = 30?) 、 
长 度 (22 = 14 mm) 对 岩石 蠕 变 破 裂 的 影响 最 具 明 
显 。 现 选 定 含 倾角 为 443" ,长 度 为 10 mm 的 单一 裂 
隙 宕 体 为 标准 岩 体 , 男 取 两 组 含 单一 裂隙 岩 体 ,只 对 
其 做 单一 变量 处 理 , 其 裂 际 尺寸 分 别 为 :A 组 , 裂 际 
倾角 为 30° ,长 度 为 10 mm; B 组 ,裂隙 倾角 为 45° ,长 
度 为 14 mm。 现 将 其 与 标准 岩 体 进行 数值 模拟 对 
比 ,结果 如 图 16 所 示 。 

通过 图 16 中 各 个 变量 参数 之 间 的 蠕 变 曲 线 对 
比 发 现 , 各 个 参数 对 减速 蠕 变 阶段 影响 不 大 , 当 步 人 
等 速 蠕 变 阶段 ,对 于 角度 变化 时 ,里 变速 率 为 1.06 x 
1078/8 ,对 于 长 度 变化 时 ,里 变速 率 为 8.24 x 10 77s, 
角度 的 变化 对 等 速 蠕 变 阶段 速率 的 影响 大 于 长 度 对 
其 的 影响 , 当 进 入 加 速 蠕 变 阶 段 之 后 ,可 以 发 现 长 度 
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对 里 变速 率 的 影响 会 大 于 角度 的 影响 ,但 长 度 变 化 
的 应 变 值 较 角度 变化 时 略 小 。 长 度 变 化 的 影响 在 于 
改变 了 岩 体 空隙 体积 , 随 着 裂 颖 长 度 的 增加 , 岩 体 空 
隙 体积 增 大 ,结构 不 均匀 性 增加 , 裂 颖 扩展 和 破裂 变 
得 更 加 严重 ;而 对 于 裂 除 角度 改变 时 ,不仅 对 岩 体 空 
隙 体积 有 影响 ,还 会 对 岩 体 整体 结构 形成 影响 , 岩 体 
随 角 度 的 改变 , 岩 体 结构 会 产生 随 角 度 而 改变 的 较 
大 剪 切 力 ,从 而 对 涯 体 破裂 产生 较 大 影响 ,因此 裂 际 
角度 的 改变 对 含 裂 际 岩 体 的 蠕 变 影响 要 更 大 。 
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Fig. 16 Comparison of creep curves of different 
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本 研究 推导 了 含 单 裂隙 岩 石 安 观 应 变 与 微 裂 纹 
长 度 、 角 度 之 间 的 关系 式 ,采用 连续 - 非 连续 软件 
GDEM 对 含 裂 陈 岩 体 蠕 变 破裂 宏 细 观 演化 规律 进行 
了 模拟 ,对 比 研究 了 不 同 裂隙 长 度 、 裂 隙 角度 下 , 含 
单 裂隙 岩 体 蠕 变 破裂 演化 规律 ,得 出 以 下 结论 。 

1) 给 出 了 恒定 载荷 下 含 单 裂隙 岩石 宏观 应 变 和 
应 变 率 与 微 裂 纹 长 度 .角度 之 间 关 系 式 ,绘制 相应 变 
化 曲线 并 进行 分 析 得 出 , 含 裂隙 涯 体 蠕 变 变形 随 裂 
隙 长度 的 增 大 而 增 大 , 随 裂 隙 角度 的 增 大 出 现 先 增 
大 再 减 小 的 趋势 ,其 中 倾角 为 30° 时 影响 最 大 。 

2) 采 用 连续 - 非 连续 软件 GDEM 对 含 裂隙 岩 体 
里 谈 破 裂 宏 细 观 演化 规律 进行 了 模拟 ,由 数值 计算 
可 以 得 出 ,岩石 蠕 变速 率 以 及 破裂 程度 随 着 裂 队 长 
度 的 增加 而 增加 , 岩 体 破裂 速度 随 着 裂隙 倾角 的 变 
大 出 现 先 增 大 后 减 小 的 情况 , 当 裂 隙 倾角 为 30 时 尤 
为 突出 ,可 见 当 裂隙 倾角 为 30° 时 岩 体 结构 容易 发 生 
破裂 。 

3) 结合 理论 分 析 和 数值 模拟 结果 , 分析 单 裂 际 
岩 体 里 变 裂 纹 扩展 演化 规律 。 裂 隙 角度 和 长 度 对 减 
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速 蠕 变 阶段 影响 较 小 ,但 角度 改变 对 等 速 蠕 变 阶段 
速率 的 影响 大 于 长 度 变 化 对 其 的 影响 , 当 进 入 加 速 
蠕 变 阶段 之 后 ,可 以 发 现 长度 对 蠕 变 速率 的 影响 会 
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大 于 角度 的 影响 ,但 是 长 度 变 化 的 应 变 值 小 于 角度 
变化 的 应 变 值 ,可 见 相 对 与 裂 际 长 度 变量 来 说 , 裂 辽 
角度 变量 对 含 裂 际 岩 体 的 里 变 影响 更 大 。 

4) 含 裂隙 岩 体 蠕 变 破 裂 裂纹 扩展 演化 规律 的 研 
究 可 为 实际 岩 体 工程 中 含 裂 际 类 岩 体 蠕 变 破裂 规律 


研究 及 蠕 变 治理 提供 理论 参考 。 
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